www.gusucode.com > 基于Matlab的MIMO通信系统仿真 含报告；司中威；了解移动通信 > 基于Matlab的MIMO通信系统仿真 含报告；司中威；了解移动通信关键技术，了解数字通信系统仿真流程，实现基本的信道编译码、调制解调等通信模块。（好评如潮，课设拿满） 学习并实现MIMO空时处理技术 学习性能分析的思路和方法/mimo/matlab for mimo 2x2/Sin_maker_ok.m

    clear,clf,close all
% 
fs = 48000;
f1 = 5000;
f2 = 11000;
N = 8*1024;
t = (0:(N/2-1))/fs;
x1 = [sin(2*pi*f1*t), sin(2*pi*f1*t)];
%x2 = [sin(2*pi*f1*t), sin(2*pi*f2*t)];

% f=linspace(-1/2,1/2,length(x1));
% X1=fftshift(fft(x1));
% subplot(211),plot(t,x1), grid
% 
% f=linspace(-1/2,1/2,length(x2));
% X2=fftshift(fft(x2));
% subplot(222),plot(f,abs(X2)), grid

Xout = [x1; x1];

Xout = Xout / (max(max(abs(Xout))));
Xout = Xout * (2^15 - 1);
% xout1=Xout(1,:);
% f=linspace(-1/2,1/2,length(xout1));
% Xout1=fftshift(fft(xout1));
% subplot(212),plot(t,xout1), grid


%Xout =0* randn(2,N);
% figure(2)
% subplot(221)
% plot(Xout(1,:))
% title('Audio Out, left')
% subplot(222)
% plot(Xout(2,:))
% title('Audio Out, right')

[Xin,no_samples_delivered,buffer_size]=DSK6713_rxtx(N,2,2,Xout,fs);

figure(2)
xin1=Xin(1,:);
xin2=Xin(2,:);
% 
% f=linspace(-1/2,1/2,length(xin1));
% Xin1=fftshift(fft(xin1));
subplot(211),plot(xin1), grid
% 
% f=linspace(-1/2,1/2,length(xin2));
% Xin2=fftshift(fft(xin2));
% subplot(212),plot(f,abs(Xin2)), grid

% % plot(Xin')
% figure(1)
% subplot(223)
% plot(Xin(1,:))
% title('Audio In, left')
% figure(1)
% subplot(224)
% specgram(Xin(2,:),256,fs)
% title('Audio In, right')
title('y(n)')
xlabel('v=fo/fs (fs=96000)')
ylabel('power')